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Abstract 


Computational analysis of two 1911 Wright brothers “ Bent End ” wooden 
propeller reproductions have been performed and compared with experimental 
test results from the Langley Full Scale Wind Tunnel. 1 The purpose of the 
analysis was to check the consistency of the experimental results and to validate 
the reliability of the tests. This report is one part of the project on the propeller 
performance research of the Wright “Bent End” propellers, intend to document 
the Wright brothers’ pioneering propeller design contributions. Two computer 
codes were used in the computational predictions. The FLO_MG Navier-Stokes 
code 2 is a CFD code based on the Navier-Stokes Equations. It is mainly used to 
compute the lift coefficient and the drag coefficient at specified angles of attack at 
different radii 3 . Those calculated data are the intermediate results of the 
computation and a part of the necessary input for the Propeller Design Analysis 
Code 4 (based on Adkins and Libeck method) , which is a propeller design code 
used to compute the propeller thrust coefficient, the propeller power coefficient 
and the propeller propulsive efficiency. 

Nomenclature 

x, y, z = physical direction 

x le , z te = the coordinate of the leading edge at the physical domain 
x te , z le = the coordinate of the trailing edge at the physical domain 


2 


z = an intermediate variable at the physical domain for the convenience of 


domain conversion 

x_ ba[ , y b „ , z_ bar = computational direction 

<J) = flow angle 
c = blade section chord 
a = angle of attack 
R = propeller tip radius 
r = radial coordinate 
= nondimensional radius = r/R 

C, = blade section lift coefficient 
C d = blade section drag coefficient 
C p = propeller power coefficient 
C T = propeller thrust coefficient 

J = propeller advance ratio 

RPM = rotational speed of the propeller 

V = forward velocity 

P = dynamic pressure 

r] = propeller propulsive efficiency 


Introduction 

The Wright brothers spent a great deal of time and resources in trying to design 
high efficient, practical wooden propellers. As we approach the 1 00th anniversary 
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of the Wright brothers’ first controlled, powered flight, it is important to recognize 
that perfecting efficient propellers was a key contribution to the Wright brothers’ 
success. Two Wright brothers bent end wooden propeller reproductions were 
built for the testing in the Langley Full Scale wind tunnel. 1 Measurements have 
demonstrated the effectiveness of their ingenious use of Wilbur Wright’s blade 
element theory, exploiting large diameter propellers, turning at low rotational 
speeds. Their “ optimized “ propeller designs utilized rearward blade sweep and 
incorporated a type of composite fiber tip covering to produce propellers with 
maximum efficiencies above 85 percent at nominal advance ratios slightly above 
one. Since the maximum efficiencies are as good as or even better than modern 
propeller designs, it is really surprising. To confirm the accuracy of the 
experiment results, computational predictions were done for comparison, using 
the FLO_MG Navier-Stokes code 2 and the Propeller Design Analysis Code 4 . 
The initial conditions are the same as the test conditions, and most of the critical 
intermediate results were checked and corrected to guarantee accuracy of the 
subsequent calculations and the final results. 

Procedure 

A. Calculations on the FLO_MG Navier-Stokes code 

The FLO_MG is a code used to acquire the blade section lift coefficient and the 
drag coefficient. 2,3 The dimensionless coordinates in the computational domain, 
the specified angle of attack, Mach numbers, and Reynolds numbers are all 


4 


necessary input for the code. The mesh plots are intermediate results, the 
accuracy of which will ensure the meaning of the later work. 

(1) Dimensionless conversion of the coordinates to computational domain 
The original geometric parameters of the propeller were obtained from the careful 
measurements at NASA Langley Research Center. They include the coordinates 
of the points at the edge of each cross section, the area of each cross section, 
and the perimeter of each section. Ten blade cross sections radii were selected 
for profile data. Their nondimensional radiuses are 0.20, 0.30, 0.4, 0.5, 0.6, 0.7, 
0.8, 0.9, 0.96 and 1.0, while the whole propeller blade radius is 51 in (4.25 ft). 
Figures 1 .01 to 1 .08 are plots of various blade cross section from | = 0.2 to 
| = 0.9 (The coordinate where § = 0.96and | = 1 is not available). Each blade 
section was represented as pairs of x and z coordinates in a data file, which is 
described by these x P z. “points”. The local chord length is computed by 

calculating all combinations of ^(X; - Xj) 2 + ( z t -z for all i and j, and 

finding the maximum length from these calculations. The relevant two points 
making the chord represented the position of the leading edge and trailing edge. 
To change the two-dimensional propeller blade cross-sections into dimensionless 
airfoils, with a chord length of unity and with the chord reference plane aligned 
with the x bar - y b> plane (z bar =0), the equations (1)-(5) were used to effect 
the conversion. 

c = V( z .c “ Z J 2 + ( x . e “O' (1) 
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<() = tan ' 1 — — — 
x te -x te 

z . =(x-x te )*tan<J) + z te 

(x k -x)/ cos - (z - z a ) *sin(|) 

X bar 

c 

(z " Z a ) * COS(j) 


■(2) 

(3) 

(4) 

(5) 


Figures 2.01 to 2.08 are the converted plots of the cross sections for £ =0.2 to £ 
=0.9 based on Equations (1)-(5). In the new plots, the dimensionless top and 
bottom surface are located at the positions according to functions of the new 
dimensionless horizontal x-coordinate (x bar)- Since in the new coordinate 
system, each section was of unit length and no longer used the airflow direction 
as the reference horizontal direction, the flow approached the section with an 
angle of attack corresponding to the actual orientation of the propeller section to 
the original coordinate system. The angles at which the air will strike each airfoil 
section had to be specified. 

(2) Reynolds number and Mach number 

Eagle Aeronautics provided the Mach number vs. Reynolds number in plot, 
shown in Figure 3.01 and Figure 3.02. From Figure 3.01 , the Mach number at 
each cross section can be gotten directly. From the Figure 3.02, the Reynolds 
number over the chord length can be calculated. Multiplying the unit Reynolds 
number by the chord length at each section results in the Reynolds number. Like 
most CFD codes, FLO_MG has stability problems when the Mach number is too 
low, as is the case for most sections. Considering the fact that compressibility 
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effects do not come into play until Mach ~ 0.3, it is safe to assume that the 
results will still be satisfactory if the Mach number is increased to computationally 
acceptable values. So Mach = 0.106 (the value for the airfoil section at | = 0.6) 
for airfoil sections at | = 0.1 to % = 0.5 was used to permit stable calculations and 
keep correct values of the respective Reynolds number. 

Besides the dimensionless geometrical coordinates, the Mach number, and the 
Reynolds number, additional information was needed to create the input file, 
such as the angle of attack, which is -4, -2, 0, 2, 4, 6, 8, and the specified 
number of grid points in computational domain, which is 321*64. 

(3) Mesh 

The mesh is the intermediate output of FLO_MG. 2 A good mesh is necessary to 
guarantee meaningful numerical results. All the input information should be 
written into the FLO_MG requested input format. There are several shapes for 
grids around two-dimensional bodies: C-grids, H-grids, and O-grids etc. C-grids 
work best when the airfoils have sharp trailing edges, and O-grids are best for 
airfoils having blunt, rounded trailing edges. Unfortunately, the grid generator of 
FLO_MG is not designed to handle airfoils like the ones in the project (O-grids), 
i.e., having a concave lower surface and a blunt trailing edge. It was designed to 
generate C-grids only. So, to get a reasonable mesh, either a different grid 
generator had to be used, or the shape of the trailing edge had to be altered. The 
latter method was chosen in this computational project to get rid of the 
generation problem: a few points around the trailing edge on both sides were 
removed to sharpen the trailing edges, while all other parameters in the input 
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files remained the same. The corresponding mesh with 321*65 grid points was 
very good. Figures 4.01 to 4.08 are the mesh plots of each cross section (from £ 
=0.2 to £ =0.9). 

(4) The blade section lift coefficient and the blade section drag coefficient 

All flow solvers, including FLO_MG, work on the principle that the relevant 
equations, usually some form of the Navier-Stokes or Euler equations, are 
iterated until a satisfactory solution is obtained. For all the project cases- 
depending on the Mach number, Reynolds number, and angle of attack -- one 
would want the solution to run long enough so that the flow around the airfoil is 
reasonably stable, i.e., it doesn’t change with time (time, of course, being 
represented by number of iterations), or the changes are periodic. One program 
of FLO_MG named exflomg is used to extract the convergence history for Ci and 
Cd. 2,3 In the first instance, one would expect a plot of, say, lift coefficient (Ci), to 
look like the ones in Figure 5.01 . Note from the first plot that the Ci fluctuates 
initially, but eventually reaches a value that remains constant regardless of the 
number of additional iterations. This solution is converged, and it represents a 
steady state for the flow (the discontinuous breaks in the convergence plot 
represent the places where the density of the grid points was increased. The 
number of iterations at which this occurs is specified in the input file for program). 
The solution in Figure 5.02 is not as steady as Figure 5.01 , but it is acceptable. In 
the second instance, the flow may become unsteady, i.e., it is time dependent. 
In this project, the prime reason for this is that the flow separated from the 
surface of the airfoil. This behavior is manifested as oscillations in the 
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convergence histories of parameters like Ci and Cd, as can be seen in Figure 
5.03 and Figure 5.04. Some times the amplitude of the oscillations remains 
constant, and sometimes it increases, as is the case in the Figure 5.04. If the 
flow is separated, then the value for the force coefficients changes with time. If 
the solution doesn’t change with time, then the output values of Ci and Cd are the 
correct values, since these are the ones corresponding to the last iteration. If the 
solution changes with time, then the output values do not necessarily represent 
the “true” values. In this case, the average of the last 150 values for Ci and C d 
found in the convergence history are used to obtain a average value steady- 
stats, which is a more appropriate representation of the true coefficient. Tables 
1.01 to 1.10 are Q and Cd vs. angle of attack at different cross sections (using 
the same values at | =0.9 for the values at § =0.96 and i =1). Figures 6.01 to 
6.10 are the relevant plots respectively (using the same values at ^ =0.9 for the 
values at l =0.96 and % =1). Figures 7.01 to 7.10 are the relevant plots 
respectively (using the same values at £ =0.9 for the values at £ =0.96 and | =1 ). 

B. Calculation on Propeller Design Analysis Code 
(1) The input 

The input should also be written into the code specified format. There are two 
input files for the Propeller Design Analysis Code. 4 The first one is a common 
input file, which is the basic input for all test cases. It includes all the Q and Cd 
gotten from FLO_MG at normalized span locations, the number of blades (which 
is 2); blade radius (which is 4.25ft); the number of specified angles of attack 
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(which is 7), and the specified angle of attack (which are -4, -2, 0, 2, 4, 6, 8), etc. 
These parameters keep constant for all test cases. The second input file is a 
special input file, which includes the propeller design RPM, design forward 
velocity, the dynamic pressure, and the relaxation parameter that can be played 
to ensure convergence. These parameters are case dependent. There are two 
primary sets of computational runs. In the first group, the Propeller RPM was 
specified as 400, with the design speed ranging from 8.69 knots to 43.45 knots. 
In the second group, the design velocity was specified at 34.76 knots, with the 
design RPM ranging from 325 to 2000. 

(2) The Output 

The thrust coefficient, Ct, power coefficient, C p , and the efficiency, r|, are the final 
output of this project. Table 2.01 supplies the calculated advance ratio J, the Ct, 
C p , and r] from the propeller code output with specified RPM equal to 400. Table 

2.02 supplies the calculated advance ratio, J, the Ct, C p , and r| from the code 
output with specified velocity equal to 34.76. Figure 8.01 shows Ct vs. J of the 
two parametric groups. Figure 8.02 shows C p vs. J of the two groups. Figure 8.03 
shows r| vs. J for the two groups of occasions. 

Analysis of the result 

Table 3.01 is a comparison of Ct between experiment and computation. Table 

3.02 is a comparison of C p between experiment and computation. Table 3.03 is a 
comparison on r\ between experiment and computation. The Figures 9.01 to 9.03 
are the relevant plots for the comparison. From the comparison, it can be seen 
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that from J=0.6 to J=1.15, the computed results complied with the experimental 
results very well. In fact, from J=0.6 to J=1.15 (the stalled advance ratio of 
experiment), the computational values are a little higher than the experimental 
ones. There are obvious differences when J is below 0.6. This is mainly because 
of the unreliability of the computational results when J is too low (The angle of 
attack in the output is far outside the specified range, which is from -4d to 8d). In 
the experiment, stall occurred at J=1.15, while in the computation the stall 
occurred at J=1 .29. In the other words, the computational ri values still increased 
from j=1 .1 5 to j=1 .29, but the experimental values decreased. In the experiment, 
max efficiency is 0.87 at j=1 .15 (in the computational result, at J=1 .15, r|=0.894); 
In the computational result, the maximum efficiency was 0.93 at j=1 .29. 
According to the empirical expectation, the stall advance ratio should not be as 
high as 1 .29 in the computation. The main reason of the computational result 
may be on the sharpen of the trailing edge in order to generate C-grids for 
FLO_MG. This alteration brought inaccuracy on the lift coefficient and drag 
coefficient, which changed the final stall position. The relaxation coefficient in the 
input file of the Propeller code may be another factor to influence the stall. The 
proper relaxation coefficient should be given other than 10 to all the cases. 
Anyway, the fault on the difference is because of some approximation in the 
computation. So it can be concluded that the experimental results from Langley 
is quite reliable, and the Wright “ Bent end” propeller does have a very high 
propulsive efficiency that is as good as or better than modern propeller. 
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Appendix: the important figures and tables 



Figure 1.01 Blade contour at ^=20% 


(Area 4.605, Perimeter 12.238, Xi, =1.742, z ta =6.615, =-0.997, =1.554, c=5.7546, «|»=61 .5779d) 


Original Plot (r/R&0.3) 



- 3-10 1 2 4 

X 


Figure 1.02 Blade contour at ^=30% 


(Area 4.1 84, Perimeter 1 3.367, x*. =2.488, z,. = 6.524 , x t . =-1 .377, z,. =1 .427 , c= 6.3967, <fr= 52.8274d) 
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Original Plot (r/R=0.4) 



- 3 - 2-101234 

X 


Figure 1.03 Blade contour at =40% 


(Area =4.330, Perimeter = 15.489, Xi. =3.441, z,. =6.483, x,. =-1.852, z t . =1.223, c= 7.4621, <fr= 44.8208d) 


Original Plot (r/R=0.5) 



-3 - 2-1 0 1 2 3 4 5 

X 


Figure 1.04 Blade contour at h, =50% 


(Area =4.647, Perimeter = 16.004, Xi. =4.39, =6.458, x t . =-2.41, z* =1.036, c= 8.6970, 4= 38.5672d) 
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Original Plot(r/R=0.8) 



-5 -4 -3 -2 -1 0 1 2 3 4 5 6 

X 


Figure 1.07 Blade contour at ^=80% 


(Area =4.170, Perimeter = 22397, Xj* =5.157, z* =5.501, x te =-4.657, =0.606, c= 10.9670, $=26.5090d) 


Original Plot(r/R=0.9) 



-6-5-4-3-2 -1 0 1 2 3 4 5 

X 


Figure 1.08 Blade contour at i; =90% 

(Area =3.350, Perimeter = 21.834, Xi. =4.28, Zi. =4.721, *. =-5.564, z«. =0.482, c= 10.7179, 4>= 23.2976) 
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Figure 2.06 Dimensionless, reoriented blade at ^ =70% 













Figure 4.02 Computational grid for ^ =30% 




Figure 4.04 Computational grid for ^=50% 





Figure 4.05 Computational grid for ^=60% 







Figure 4.08 Computational grid for 1= =90% 
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Figure 5.03 Representative C, convergence history for an unstable convergence, 

showing a growing oscillation 



Figure 5.04 Representative C d convergence history for an unstable convergence 
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Figure 6.01 Lift coefficient VS. angle of attack at H, =20% 



Figure 6.02 Lift coefficient VS. angle of attack at H, =30% 
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Figure 6.07 Uft coefficient VS. angle of attack at =80% 


Figure 6.08 Lift 
coefficient VS. angle 
of attack at S;=90% 
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Figure 6.09 Lift 
coefficient VS. angle 
of attack at £; =96% 
(Same as H, =90%) 


Figure 6.10 Lift 
coefficient VS. angle 
of attack at £*=100% 
(Same as | =90%) 
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Figure 7.03 Drag coefficient VS. angle of attack at $5 =40% 



Figure 7.04 Drag coefficient VS. angle of attack at =50% 
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Figure 7.06 Drag coefficient V 
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Figure 7.07 Drag coefficient VS. angle of attack at =80% 


drag coefficient VS. angle of attack 
(r/R=0.9) 



- 6-4 -2 02468 10 


Figure 7.08 Drag coefficient VS. angle of attack at ^ =90% 
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drag coefficient VS. angle of 
attack (r/R=0.96) 



- 6-4 -2 02468 10 


Figure 7.09 Drag coefficient VS. angle of attack at ^ =96% 
(Same as I; =90%) 


drag coefficient VS. angle of 
attack (r/R=1) 



- 6 - 4-2 02 4 68 10 


Figure 7.10 Drag coefficient VS. angle of attack at ^ =100% 
(Same as ^ =90%) 
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Figure 8.01 Comparison of fixed rotational speed and fixed flow-speed CFD 
predictions of thrust coefficient VS. advance ratio 



Figure 8.02 Comparison of fixed rotational speed and fixed flow-speed CFD 
predictions of power coefficient VS. advance ratio 
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The efficiency VS advance ratio 



Specified RPM=400 — • — Specified Vek>city=34.76 knots 


Figure 8.03 Comparison of fixed rotational speed and fixed flow-speed CFD 
predictions of efficiency VS. advance ratio 


CTVSJ 



Experimental Values Computational Values 



Figure 9.01 Comparison between numerical predictions and experimental 
measurements of thrust coefficient VS. advance ratio 
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Figure 9.02 Comparison between numerical predictions and experimental 
measurements of power coefficient VS. advance ratio 
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Figure 9.03 Comparison between numerical predictions and experimental 
measurements of efficiency VS. advance ratio 
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Table 1.01 Lift coefficient and Drag coefficient VS. angle of attack at ^=20% 


( Mach=0.106, Rey= 206000) 


a 

c, 

c d 

HSBES 


2.25E-02 


3.34E-01 


O.OOE+OO 

5.63E-01 

2.20E-02 

2.00E+00 

7.33E-01 

2.32E-02 

4.00E+00 


2E2SE 

6.00E+00 

1.1 IE-i-00 

2.70E-02 

8.00E+00 

1.29E+00 



Table 1.02 Lift coefficient and Drag coefficient VS. angle of attack at ^=30% 

(Mach=0.106, Rey= 267000) 


a 

Ci__ . 

c d 

■4.00E+00 

2.55E-01 

2.16E-02 

fr*«i»3=n»n; 

4.64E-01 


0.00E+00 

6.65E-01 

2.01 E-02 

2.00E+00 



4.00E+00 

1.05E+00 


6.00E+00 

1.24E+00 


8.00E+00 

1.41E+00 

2.81 E-02 


Table 1.03 Lift coefficient and Drag coefficient VS. angle of attack at %=40% 

(Mach=0.106, Rey=354000) 
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Table 1.04 Lift coefficient and Drag coefficient VS. angle of attack at I; =50% 

(Mach=0.106, Rey= 478000) 



.00E+ 


0.00E+ 


2.00E+ 


4.00E+ 


6.00E+ 


8.00E+ 


2.79E-01 

4.07E-0 

5.35E-01 

2.02E-0 


9.47E-01 1.86E-0 


1.14E+00 1.99E-0 


1.32E+0Q 2.18E-0 


1.49E+0Q 2.49E-0 



Table 1.05 Lift coefficient and Drag coefficient VS. angle of attack at£* =60% 

(Mach= 0.106, Rey= 625500) 



(Mach=0.121, Rey= 794000) 
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a 

Ci 

c d 

-4.00E+00 

1.77E-01 

2.91 E-02 


4.12E-01 

1.60E-02 

0.00E+00 

6.24E-01 

1.53E-02 

2.00E+00 



4.00E+00 

1.03E+00 

1.68E-02 

6.00E+00 

1.21E+00 

2.02E-02 

8.00E+00 

1.32E+00 

3.00E-02 


Table 1.07 Lift coefficient and Drag coefficient VS. angle of attack at ^=80% 

(Mach= 0.134, Rey= 868000) 


a 

c, 

c d 

■4.00E+00 

1.25E-01 

2.85E-02 



■ir^»=g.>j 

O.OOE+OO 

5.82E-01 

1.47E-02 

2.00E+00 

7.93E-01 

sbhe 

4.00E+00 

9.92E-01 

1.64E-02 

6.00E+00 

1.16E+00 

2.34E-02 

8.00E+00 

1.26E+00 

3.83E-02 


Table 1.08 Uft coefficient and Drag coefficient VS. angle of attack at £• =90% 

(Mach= 0.148, Rey= 938000) 


a 

c, 

Cd 

-4.00E+00 

1.25E-01 

2.85E-02 


3.65E-01 

1.61 E-02 

0.00E+00 



2.00E+00 

7.93E-01 

1.51 E-02 

4.00E+00 

9.92E-01 

1.64E-02 

BE5BES 

1.16E+00 

2.34E-02 

8.00E+00 

1.26E+00 

3.83E-02 


Table 1.09 Uft coefficient and Drag coefficient VS. angle of attack at ^=96% 


(Ci, and C d are same to £ =0.9) 
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a 

c, 

c d 

■4.00E+00 



BUSHES 


rrrrrr^ 

0.00E+00 

5.82E-01 

1.47E-02 

2.00E+00 

7.93E-01 

1.51E-02 

4.00E+00 

9.92E-01 

1.64E-02 

6.00E+00 

1.16E+00 

2.34E-02 

8.00E+00 

1.26E+00 

3.83E-02 


Table 1.10 Lift coefficient and Drag coefficient VS. angle of attack at£; =100% 

(Ci, and C d are same to £ =0.9) 



Dynamic 

Pressure 

J 

? 

Ct 

Cp 


0.255941 

0.258753 



tiEEm 

17.378 

1.023763 

EEES5 


K iWi L 

0.174406 

21.7225 

1.599629 

umw 



0.176556 


2.303466 

usnm 

ncmrr 

0.168139 

msm 

WHERE 

3.135273 

C« 

0.813726 

0.145245 

0.16165 

34.756 

4.095051 

1.035012 

0.860441 

0.120459 

0.144898 

36.92825 

4.622928 

1.0997 


0.107695 

0.134429 

39.1005 

5.182798 

BEES 

1 A’AVJFmh 


WWW, 

41.27275 

5.774661 

1.229077 

0.9182 

0.080768 


bebe 

mmri 

gra-gy/s- 

lli < 

r mm 

0.066127 



Table 2.01 CFD prediction of thrust coefficient, power coefficient and efficiency at 
fixed rotational speed (Specified RPM=400) 
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RPM 


Dynamic 

Pressure 


2000 


4.095051 


0.207100 0.3199 


0.243558 0.157681 


1500 


4.095051 


0.276141 


0.400486 


0.236301 


0.162939 


1000 


4.095051 


0.414212 0.53247] 


0.220589 0.171597] 


700 


4.095051 


0.591731 


0.659026 


0.197824 


0.177624 


600 


4.095051 


0.690353 0.71458 


0.183564 


0.17734 


500 


4.095051 


0.828424 0.780640 


0.160659 0.170492 


400 


4.095051 


1.035529 0.859551 


0.121736 


0.146650 


375 


4.095051 


1.104565 0.8814A 


0.108042 0.135387] 


350 


4.095051 


1.183462 0.904620 0.091789 0.120082 


325 


4.095051 


1.274490 0.929370 0.071972 0.098699 


Table 2.02 CFD prediction of thrust coefficient, power coefficient and efficiency at 
fixed flow speed (Specified Veloclty=34.76 knots) 


J 



0.2 

0.22 

0.244 

0.4 

0.198 

0.221 

0.6 

0.185 

6.192 


0.158 

0.162 

i 

0.122 

0.126 

i.i 

0.102 

0.109 

1.15 

0.092 

0.095 

1.2 


0.085 

1.3 

0.05 

0.064 


Table 3.01 Comparison between numerical predictions and experimental 
measurements of thrust coefficient VS. advance ratio 
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measurements of power coefficient VS. advance ratio 



measurements of efficiency VS. advance ratio 
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